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24 I. INTRODUCTION 



25 The r-process likely happens in supernovae [1] , in neutron star mergers [2] , or some other 

26 high temperature (T > 1 GK) and high neutron flux (0 > 10^*^ cm^^s^^) environment. At 

27 present, the case for an indisputable r-process site has not been made, and recent arguments 

28 point out the necessity of multiple sites [3]. Supernovae have long been cited as potential 

29 r-process factories because they produce the necessary explosive conditions and they occur 

30 frequently enough to produce substantial abundances. The production of heavy nuclides 

31 arising from explosive nucleosynthesis at a Type II supernova site is linked to the rate of 

32 ^Be production [4]. 

33 The behavior of a core-collapse supernova has been described, for example, by Woosley 

34 and Janka [5]. A star of 8 to 25 M© passes through the stages of hydrogen, helium, carbon, 

35 neon, oxygen, and silicon burning at its core, continuously growing hotter and more dense. 

36 It eventually forms an iron-group core of Chandrasekhar mass (!«1.4 M©). Since the core 

37 has no other source of thermonuclear energy to support the pressure, it collapses. Energy 

38 is released in the form of neutrino radiation, while the collapse is accelerated by electron 

39 captures and photodisintegrations. An Earth-sized iron-core collapses with a velocity of 

40 ~ c/4 into a single, neutron-rich nucleus about 30 km in diameter. The collapse generates 

41 a rebounding shock wave, which ultimately stalls as it attempts to push through the in- 

42 falling matter. The proto-neutron star briefly continues accreting matter, while radiating 

43 ~ 10^^ J in the form of neutrinos, accounting for nearly 10% of its rest mass. At this stage, 

44 convection, rotation, and magnetism likely contribute to the dynamics of the subsequent 

45 explosion. However, self-consistent, three-dimensional models that predict explosion remain 

46 elusive [6]. 

47 A. The r-process 

48 The r-process produces about half of the nuclides heavier than iron [5]. It requires a hot, 

49 neutron-dense environment where neutron captures occur so rapidly that the nucleosynthesis 

50 path is pushed far out to the neutron-rich side beyond the stability valley. After cessation 

51 of the neutron flux, the short-lived nuclei /3-decay to stable species. 

52 Preceding the r-process is the a-process [7] , which is driven by charged-particle reactions 
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53 and takes place when the shock has cooled from 5 GK to 3 GK over a time period of a few 

54 seconds. In this stage of nucleosynthesis a reaction path is needed to bridge the stability 

55 gaps at A = 5 and A = 8; the most efficient path is a + a ^Be and ^Be (n, 7) ^Be, followed 

56 by^Be + a ^^C + n [8]. As cooling continues this reaction sequence largely establishes the 

57 neutron-to-seed- nucleus ratio to which the subsequent r-process is very sensitive [9]. Too 

58 few seed nuclei will under-produce r-process nuclei, while too many seed nuclei produced in 

59 the a-process will starve the r-process environment of neutrons. Since r-process abundance 

60 predictions in certain stellar models are extremely sensitive to the a (an, 7)^Be rate [4, 10, 

61 11], establishing a precise rate for the formation of ^Be via the a(an,7)^Be reaction is 

62 required for accurately modeling nucleosynthesis in supernovae. 

63 The ^Be (7, n) ^Be reaction may be used to deduce the ^Be (n, 7) ^Be cross-section by 

64 applying the reciprocity theorem. Previous ^Be photodisintegration studies are numerous 

65 [12-26] but reveal relatively large cross-section uncertainties in the astrophysically important 

66 region near the neutron emission threshold. New (7,^,) cross-section measurements with 

67 improved accuracy are now possible using intense photon beams with small energy spreads 

68 and neutron detectors with large solid angle coverage and high efficiencies. 

69 In the following. Sect. II describes a new measurement of ^Be (7, n) ^Be and provides de- 

70 tails of the data analysis used to obtain cross sections. Section III describes the methods used 

71 to extract resonance parameters from the new data using proper energy dependences for par- 

72 tial widths of ^Be excited states. The methods employed to calculate the ^Be (n, 7) ^Be cross 

73 section and the corresponding a{an,'-f)^Be reaction rates are described in Sect. IV, along 

74 with a comparison of the determined rates with those from earlier studies. Finally, Sect. V 

75 presents a summary of the findings presented in this paper. 

76 II. EXPERIMENT 

77 CoUimated, near-monoenergetic photon beams of 1.5 MeV < < 5.2 MeV were incident 

78 on a thick ^Be target located within the central bore of the neutron detector. The abso- 

79 lute number of neutrons from the (7,^,) reaction was determined using a moderated ^He 

80 proportional counter with a high efficiency (--60%) for detecting low-energy neutrons [27]. 

81 The absolute incident photon flux was measured using a large Nal(Tl) detector. Photon 

82 beam energy resolution was determined with a high purity germanium (HPGe) detector. A 
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FIG. 1. Schematic diagram of the experimental setup for the ^Be (7, n) '^Be reaction measurements. 
After collimation, the photon beam passes through scintillation paddles (not shown) and into 
the target room. The photon beam then passes through the following elements: (a) "clean-up" 
collimator wall; (b) the chosen target located near the longitudinal center of the neutron counter; 
(c) lead attenuators located between lead collimator walls; (d) Nal(Tl) detector; (e) HPGe detector. 

83 schematic diagram of the experimental setup is shown in Fig. 1. 

84 A. Experimental setup 

85 Intense, collimated photon beams (0 r; 3 x 10^ 7/s) are routinely produced at the Tri- 

86 angle Universities Nuclear Laboratory's High Intensity Gamma Source (HI7S) by inverse- 

87 Compton scattering of free-electron-laser photons from electron bunches circulating in a 

88 storage ring [28]. For the present experiment, a circularly-polarized beam was used and some 

89 flux was sacrificed to attain the high photon energy resolution needed to map the detailed 

90 behavior of the cross section at the three-body (1573 keV) and two-body (1665 keV) thresh- 

91 olds shown in Fig. 2. Present data were taken using beam intensities of 10^ < (f) < 10^ 7/s and 

92 energy spreads of < 1%. 

93 The photon beam was defined by a 12-mm diameter, 30.5-cm thick lead collimator. It 

94 then passed successively through three thin scintillating paddles that acted as a relative 

95 photon flux monitor, a 2.54-cm diameter hole in a lead "clean-up" collimator (CC), and 

96 -'1.5 m of air before reaching a second CC (shown as (a) in Fig. 1) placed directly in front 

97 of the neutron detector. 

98 At each energy, photons impinged successively on one of three thick, 19-mm diameter 

99 cylindrical targets (^Be, D2O, graphite) described in Table I, or air, as they passed through 
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FIG. 2. (Color Online) A level scheme relating the mass energies of a + a + n, ®Be + n, and 
^Be. Thresholds for three-body and two-body breakup of ^Be are shown to occur at incident 7-ray 
energies of 1573 keV and 1665 keV, respectively. The latter is also the neutron separation 
energy Sn for ^Be. Energies shown for ^Be excited states are from the present work. In the rate 
calculation (described in Sect. IV), E is the center-of-mass energy of the two a-particles. The 
parameter E' is the energy of the ®Be nucleus and the neutron with respect to E. In this scheme, 
formation of ^Be at ii^ = £" = is very unlikely, but not prohibited because the ground state of 
^Be has finite width. 



TABLE I. Targets used in the present experiment and their physical properties. 



Material 


Length Density 
(cm) (g cm"'^) 


Mass 
(g mor^ 


Thickness 
) (nuclei cm"^) 


^Be 


2.54 


1.848 


9.012 


3.14x10^^ 


Graphite 


2.54 


1.700 


12.01 


2.89x10^3 


D2O 


7.59 


1.106 


20.04 


5.05x10^3 
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FIG. 3. (Color Online) A HPGe spectrum for 2470-keV photon beam with a resolution of AE/E = 
1%. Spectra obtained with ^''Co, ^"^K, and "^^^Tl calibration sources are overlaid. Images showing 
confirmation of target alignment are inset, (a) An unattenuated beam profile. Flux was nearly 
constant across the 12-mm diameter and decreased rapidly at the edges, (b) Contrast from a 4-mm 
diameter lead alignment pellet confirmed the axial placement of the target. 

100 the neutron detector. The heavy water target was bombarded under the same experimental 

101 conditions as the ^Be target, allowing for normalization of the ^Be (7,n) ^Be measurements 

102 to the well-known {'~f,n) cross section [29] and for calibration of the neutron detector 

103 efficiency. The graphite target was used to determine beam-induced backgrounds in the 

104 neutron detector. 

105 To increase the efficiency of data collection, targets were remotely rotated into the beam 

106 using a four-position Geneva mechanism, which also assured reproducible alignment of the 

107 axis of each cylindrical target sample with the beam axis. Axial alignment of the target was 

108 confirmed using an alignment pellet and a photon beam imaging system [30] (see Fig. 3). 

109 Downstream of the target, lead of various thicknesses could be inserted to attenuate the 
no beam by up to a factor of 100. Further downstream, the remaining photons were incident on 

111 either a Nal(Tl) detector or a HPGe detector, depending on whether photon flux or energy 

112 was being measured. The lead attenuators facilitated simultaneous high neutron counting 

113 rates and negligible Nal(Tl) signal pile-up. Data acquisition dead-times were assured to be 

114 small. The counts measured in the Nal(Tl) detector were corrected for detection efficiency 

115 as well as for attenuation through the lead and targets. 
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FIG. 4. (Color Online) Effective neutron detector efficiency versus photon energy for 
^Be (7, n) ^Be with error bars representing the 4.6% systematic uncertainty. Above E^y - 2431 keV, 
average neutron energy, and thus detection efficiency, was determined using the experimental IjO 
ratio. The (red) curve is modeled using Eq. 1 with the branching ratios listed in Table II. Discrep- 
ancies in the 3 to 4 MeV range are not well understood, but do not affect the present calculation 
for the astrophysical a(an,7)^Be reaction rate discussed in Sect. IV. 

B. Detector calibration 



117 Absolute measurements of the number of photons on target and the number of emitted 

118 neutrons from the reaction were needed to determine the total cross section of ^Be (7, n) ^Be. 

119 Thus, it was essential to determine the absolute energy- dependent detection efficiencies of 

120 the neutron detector and the large Nal(Tl) detector. The active neutron detection elements 

121 were 18 tubular proportional counters, each containing ^He at 6.1 x 10^ Pa. The tubes, 

122 embedded in a cylindrical polyethylene body that served as a neutron moderator, were 

123 arranged in concentric inner (/) and outer (O) rings of nine equally spaced detectors each. 

124 The energy- dependent efficiency of the neutron counter was determined in an extensive 

125 study [27]. The ratio of counts in the inner and outer rings {IjO ratio) provided a coarse 

126 estimate of the average neutron energy. 

127 The total efficiency of the large Nal(Tl) detector was found to be nearly constant (98.3 ± 

128 1.7%) over the experimental energy range using the Monte-Carlo codes geant4 and mcnpx. 

129 The results were consistent with data obtained using the ^^F(p, 0:7) reaction, taken using a 

130 mini-tandem accelerator [31]. This measurement provided a determination of the absolute 
132 detection efficiency for 6.13 MeV photons [32, 33]. 



133 While the efficiency of the large Nal(Tl) detector was constant for photons within the 

134 experimental energy range, the effective efficiency of the neutron detector for photoneutrons 

135 from ^Be, e^siE^), varied and is shown in Fig. 4. As can be seen in Fig. 2, the average 

136 energy for neutrons which decay to the ground state of ^Be is Ej-Sn, and for these neutrons 

137 the detection efficiency e„ may be described using a sixth-order polynomial in En- Above 

138 Ej = 2431 keV, multiple decay channels exist and the average neutron energy is no longer 

139 proportional to the photon energy. Instead, it becomes nonlinear and strongly dependent 

140 on the photon energy, which causes tcs to deviate from the polynomial that describes e„. To 

141 account for this, the effective efficiency was described using 

.„(^£„)(E^i^^). (1) 

\ j ^tot I 

142 where for level j, cTj is the contribution from state j to the total cross section atot-, and fij is a 

143 branching ratio for state j to the ground state of *Be. This form assumes that a newly opened 

144 neutron branch will decay with a small neutron energy 5 En- The detection efficiency for the 

145 fraction of the total cross-section decaying to the ground state is described by the simple 

146 polynomial e„ {E^ - Sn), while the efficiency of new branches is e„ (6En) ~ emax- This model 

147 is compared in Fig. 4 to a point-by-point determination of the neutron detector efficiency 

148 constructed using the I/O ratio. The branching ratio for the 5/2" state at 2431 keV was 

149 taken from Ref. [34]. Other branching ratios were chosen to make the model congruent 

150 with the point-by-point analysis. The contribution of states other than the 1/2+ state to 

151 the a {an,'j)^'Be rate will be shown in Sect. IV to be nearly independent of the choice of 

152 branching ratios. 



153 C. Data analysis 

154 We performed two analyses of the yield data: (i) assuming a monoenergetic photon beam; 

155 and (ii) assuming a photon beam with a finite energy width, requiring deconvolution of the 

156 photon beam energy profile to interpret the neutron yield accurately. 

157 Under the assumption of a monoenergetic beam, the cross section may be written as 

N,-(NTlA)-en' 
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158 where is the number of detected neutrons, is the number of incident photons, Nt/A is 

159 the effective number of target nuclei per unit area, and e„ is the neutron detector efficiency. 

160 The quantity Nrp/A was determined by comparing neutron yields from target-in and target- 

161 out runs. Thick targets required a photon energy- dependent correction of the form 



162 where /i is a material specific attenuation coefficient and t is the thickness of the target. 

163 This correction accounted for the reduction in the number of incident photons caused by 

164 interactions within the target volume. 

165 Total ^Be(7,n)^Be cross-section uncertainties for a in Eq. 2 were found to be 3.2% 

166 (statistical) and 4.6% (systematic). The largest contributions to the uncertainties came 

167 from the absolute efficiencies of the neutron and Nal(Tl) detectors and the exact photon 

168 beam flux loss associated with lead attenuators. Numerical values for cross sections and a 

169 detailed analysis of the associated experimental uncertainties are available in Ref. [32]. 

170 The second analytic approach treated the reality that the photon beam was not truly 

171 monoenergetic. Such a treatment is especially important near threshold, where the cross 

172 section changed significantly within the energy spread of the beam (see Fig. 5). The 

173 experimental yield Y may be defined as 



174 where / is an energy-dependent function describing the energy distribution of the pho- 

175 ton beam and at is a trial cross section. 

176 To determine /, the detector response function was deconvolved from the HPGe spectrum 

177 at each beam energy. The resulting spectra were then normalized such that 



178 where A^^ was determined using the Nal(Tl) detector. The trial cross section at was assumed 

179 to be the sum of six Breit-Wigner equations (BWEs), each of which had three free parameters 

180 corresponding to the resonance energy Er, the neutron partial width F^, and the transition 

181 strength B(E1) or B(M1) of an excited state in ^Be. 

182 Histograms were constructed for the other components of Y{a, e„, Nt/A) as a function of 

183 photon energy. Each bin of the histogram then represented the respective component of the 



Y 



(4) 




(5) 
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184 yield over dE, the width of the bin. Thus, the calculated yield y*was given by 

i 

185 where the width of the i^^ bin was -1.6 keV. In this way the yield-weighted effective energy, 

186 E* was defined as, 

YE,J,-areMIA), 

i 

187 The trial cross section was then iteratively adjusted (over ~ 8 steps) until the global devia- 

188 tion between the calculated yield and the experimental yield was minimized. This process 

189 resulted in a deconvoluted cross section shown in Fig. 5 along with that deduced from 

190 the monoenergetic beam analysis. Figure 6 shows a histogram of the deviations between 

191 experimental yields and yields calculated from the deconvoluted cross section for 52 data 

192 points. The relative error is within ±5% at the 68% confidence interval and ±10% at the 
198 95% confidence interval. 



196 D. The Be (7, n) Be cross section 

197 From the two-body threshold energy of Eth = Sn = 1665 keV to 1900 keV, the present 

198 cross sections are larger than most of the previously reported data. Notice in Fig. 5 that 

199 deconvolution changes the shape of the cross section at the threshold, transferring strength 

200 from below to above Eth- The need for deconvolution could well explain some disagreement 

201 in this region with yields measured in earlier experiments using photon beams having larger 

202 energy spreads than those used in the present experiment [22, 25]. 

203 The narrow 5/2" state at 2431 keV was far better resolved in the present experiment 

204 than in previous works. The present experimental yield peaks more than a factor of 3 above 

205 the peak measured by Ref. [25] and nearly a factor of 2 above that of Ref. [22]. The present 

206 data, over the broad 1/2" and 5/2+ resonances near 3 MeV, are in fair agreement with the 

207 data of Refs. [18, 22] but not in agreement with the results of Ref. [25] which were obtained 

208 with a large photon beam energy spread. 

209 At energies above the broad peak at 3 MeV the present data agree with Refs. [18, 25], 

210 but not with Ref. [22]. A 3/2+ state at 4.7 MeV and a 3/2" state at 5.6 MeV are the next 
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FIG. 5. (Color Online) A plot of monoenergetic beam cross sections (black dots) and deconvo- 
luted cross sections (open red squares). Error bars are shown when uncertainties are larger than 
the points. The uncertainty shown for the monoenergetic beam data is purely statistical, while 
that shown for the deconvoluted data also includes statistical and systematic flux uncertainties 
associated with energy binning near threshold, (a) The full experimental range. The resonance 
at 2431 keV was too narrow for accurate deconvolution. Several points with high uncertainty, re- 
sulting from low-statistics runs, are apparent at ~2950 keV; these data, taken at the same energy, 
were combined for the purpose of deconvolution. (b) An expanded view of the boxed region near 
threshold better shows the differences between the two methods. The dotted (dashed) vertical line 
denotes the three-body (two-body) threshold. The sample photon beam profile shown (blue) is 
peaked at 1625 keV with an energy spread which was typical of this experiment. 



211 known excited states in ^Be [35]. These are broad states [36] which decay more strongly 

212 through the 2+ excited state in ^Be. 
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FIG. 6. (Color Online) Histogram of the relative deviation between Y and Y* . The Gaussian fit to 
the error is centered just below zero. The ±la interval is darkly shaded (teal) and the ±2a interval 
is shaded lightly (blue). 
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FIG. 7. (Color Online) (a) Data for the total ^Be(7,n)^Be cross section collected using several 
different 7-ray sources, including virtual photons from inelastic electron scattering [26] and real 
photons from both inverse Compton scattering [25] and natural radioisotopes [12, 14, 19]. (b) An 
expanded view of the boxed region better shows the different evaluations near the threshold. 
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213 III. THRESHOLD BEHAVIOR AND ENERGY DEPENDENCE 



214 Previous measurements of the direct ^Be+7 a + a + n three-body reaction in the energy 

215 range 1570 keV < < 1670 keV have yielded only an upper limit to the total cross section of 

216 93 nb [37], and the present experiment was not sufficiently sensitive to improve this. Above 

217 the two-body threshold at 1665 keV, the total ^Be(7,n) cross section rises rapidly to exceed 

218 1 mb. This feature is most logically attributed to the newly-opened decay channel, and 

219 thus excitation of the broad 1/2+ state located immediately above the two-body threshold 

220 is most frequently followed by decay to ^Be^s + n. 



221 A. Behavior near threshold 



222 Any resonance close to threshold experiences a distortion of its normal Lorentzian cross- 

223 section shape. When a level of spin J is isolated from other levels of the same spin and 

224 parity, a one-level R-matrix approximation may be used to describe the contribution of this 

225 level to the (7, n) cross section. For (7, n) reactions, this takes the form of the BWE for 

226 an isolated resonance [38]: 

a,,„^i^,; - ^2 2(2/ + 1) (^^ - E^y + ' ^ ^ 



227 with I the spin of the target nucleus and k'^ given by 

'E, 



7 

2 



228 The neutron partial width is generally written as [39] 

r„ = 272P,, (10) 

229 where 7^ is the reduced width and Pe is the penetration factor. The reduced width incor- 

230 porates the unknown parts of the nuclear interior while Pg is completely determined by the 

231 conditions outside the nucleus and may be written as 



232 where R is the channel radius, k is the wave number, and i is the neutron orbital angular 

233 momentum [38]. The channel radius R is defined as [38] 

R = ro{Af + Af), (12) 
13 



234 with At and Ap the mass numbers of the target and projectile, respectively, and tq = 1.44 fm. 

235 For neutrons, the Coulomb wave functions, Fi and Gg, are related to spherical Bessel (je) 

236 and Neumann (ne) functions by Fi = {kr)j^[kr) and Gn = (kr)ni(kr). In the cases of 

237 (7,ri) and (n, 7) reactions, the penetration factors can be written analytically, and for < 

238 i<3,Pi is [39] 

Po - kR = ./Wn. (13a) 



p - (^^»)^^^ nsci 

P (^^»)^^^ nsd) 

' 225 + 45e^„ + 6(e^02 + (^^^)3- ^ > 

239 In these expressions R is again the channel radius defined in Eq. 12; the neutron energy is 

240 related to the photon energy by = - Sn', and we define C, = 2fiR'^h''^, where fi is the 

241 reduced mass of ^Be + n. 

242 One must include Pi energy-dependence for the 1/2+ threshold resonance to obtain a 

243 good fit to the ^Be (7,n) ^Be cross-section data. However, previous works have not included 

244 energy dependence in the tails of the broad, higher-lying states in ^Be. The result of this 

245 incomplete treatment has been to infiate the previously-deduced off-resonance contributions 

246 to the a {an, 7) ^Be rate by as much as a factor of five. The states in ^Be excited by an 

247 L = 1 photon are coupled to the ground state of ^Be through emission of a neutron with 

248 a specific orbital angular momentum i determined by the spin and parity of the excited 

249 state. For the excited states considered in this paper, those with J'^ = 1/2+ , 1/2", 3/2+, 

250 3/2", 5/2+, and 5/2", decay to the ground state of ^Be through emission of a neutron with 

251 i = 0, 1, 2, 1, 2, and 3, respectively. The value of i determines the form of Pi (Eqs. 13a - 

252 1 3d), which in turn determines the energy dependence of the neutron partial width (defined 

253 in Eq. 10) for the excited state in ^Be, and ultimately the behavior of its cross section 

254 near the threshold. Figure 8 displays the relative contributions of the 1/2+ , 1/2", and 5/2+ 

255 resonances to the total cross section. The latter two are shown with and without the proper 

256 Pi energy dependence. 
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FIG. 8. (Color Online) Example of the differences between energy-independent and energy- 
dependent calculations of the ^Be(7,n)*Be cross sections at threshold. The dashed (black), dot- 
dashed (red) and dotted (blue) lines show the 1/2+ , 1/2", and 5/2+ states, respectively. Energy- 
dependent cross-section contributions from the 1/2" and 5/2+ states are identified and contrasted 
with energy-independent determinations of these states; at threshold, the latter are orders of mag- 
nitude too large. 

The energy-dependent 7-ray partial widths may be cast in terms of reduced transition 
probabilities [39]. For El and Ml transitions one finds 

r^(^l) = ^a{hc)-'E'^BiEl) i, (14) 

and 

(15) 



F^(Ml) = ^a{2MpC^)-^E^^B{Ml) i, 

where a is the fine structure constant. Note that the strength of a transition from the ground 
state to an excited state (Bf) is related to the strength of transition from that excited state 
to the ground state (Bj) by 

9 / 4- 1 

(16) 



2Jo + 1 
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TABLE II. Resonance parameters and neutron branching ratios from the present work. The latter 



were determined from the neutron detector efficiency analysis discussed in Sect. II B 



J"" 




Eb (MeV) 


El- (e^fm^) 
Ml -> ialj) 


(eV) 


Tn (keV) 


Pi (%) 


1/2+ 


El 


1.731 ±0.002 


0.136 ±0.002 


0.738 ±0.002 


213 ±6 


100 


5/2" 


Ml 


2.431 ±0.004 


0.587 ±0.027 


0.098 ±0.004 




6^ 


1/2- 


Ml 


2.880 ±0.016 


6.5 ±0.7 


1.8 ±0.2 


393 ± 18 


100 


5/2+ 


El 


3.008 ±0.004 


0.016 ±0.002 


0.45 ±0.07 


163 ± 15 


70 


3/2+ 


El 


4.704^ 


0.068 ±0.007 


7.8 ±0.4 


1541 ± 115 


38 


3/2- 


Ml 


5.59'^ 


7.8 ±2.1 


15.7±4.2 


941 ± 164 


38 



This value could not be obtained using the present data. 
^ This value was fixed in accordance with Ref. [35]. 
This was fixed in accordance with Ref. [34]. 



263 Transitions from the ground state of ^Be, Jq = 3/2", to an excited state, = 1/2, 3/2, or 

264 5/2, yield ratios of 0.5, 1, and 1.5, respectively, in Eq. 16. Since B i are the strengths of 

265 transitions used in the calculation of the a {an, 7) ^Be rate, they will be used exclusively in 

266 the following discussion. For each resonance, the reduced width 7^, the transition strength 

267 B(El)i or B(Ml)i, and the resonance energy Er are determined by fitting the data. Table II 

268 displays the parameters determined for each resonance. The sub-eV width of the narrow 5/2" 

269 state precluded experimental determination of the associated r„, and therefore the width 

270 reported in Ref. [35] was adopted for the present rate calculation. Known contributions 

271 from B(E2)| for the negative parity states were negligible compared to B(Ml)i [35]. 

272 B. Narrow resonance treatment 

273 In the case of narrow resonances, may be deduced by integrating the cross section. 

274 For this to be valid: (a) the resonance must be sufficiently isolated from other resonances, 

275 (b) the neutron and 7-ray partial widths must be small enough to be considered energy- 

276 independent, and (c) the neutron partial width must be much larger than the 7-ray partial 

277 width such that F„ !« F. With these three conditions satisfied, Eq. 8 may be integrated for 
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TABLE III. Resonance parameters and reduced transition probabilities of the 1/2^ state of 
^Be from virtual and real photon experiments. Refs. [20, 24, 40] include reanalyses of data 
originally analyzed by Refs. [19, 21, 25], respectively. 



Reaction 


Rcf. 


Er (MeV) 


Tn (keV) 


(eV) 


B(El)i (e^fm^) 


(e,e') 


Clerc et a/. [16] 


1.78 


150 ±50 


0.3 


0.050 ±0.020 


(e,e') 


Kuechler et al. [21] 


1.684 


217 ± 10 


0.27 


0.054 


(e,e') 


Click et a/. [41] 


1.68 


200 ± 20 


0.34 


0.068 


(e,e') 


Barker et a/. [24] 


1.732 


270 


0.75 


0.137 


(e,e') 


Burda et a/. [26] 


1.748 


274 ±8 


0.302 ±0.045 


0.054 


{l,ri) 


Barker et a/. [20] 


1.733 


227 ±50 


0.577 


0.106 ±0.018 




Angulo et a/. [42] 


1.731 


227 ± 15 


0.51 ±0.10 


0.094 ±0.020 




Utsunomiya et a/. [25] 


1.748 


283 ± 42 


0.598 


0.107 ±0.007 




Sumiyoshi et aZ.[40] 


1.735 


225 ± 12 


0.568 


0.104 ±0.002 




Present 


1.731 ±0.002 


213 ±6 


0.738 ±0.002 


0.136 ±0.002 



278 the 5/2 resonance, yielding 

/ a,,XE,)dE, = l(^^jT,. (17) 

279 The value obtained from this analysis is shown in Table II. 

280 C. The 1/2+ threshold resonance 

281 As will be shown in Sect. IV, the 1/2+ threshold resonance is the largest contributor to 

282 the a (an, 7) ^Be reaction rate. Table III contains resonance parameters for this state from 
288 several works. Notice that all but one of the evaluated virtual photon (e, e') data produced 

285 reduced transition strengths and 7-ray partial widths that are about half of the value of their 

286 real photon (7, n) counterparts. The anomalous parameters within the (e, e') subset of the 

287 data [24] result from a reanalysis of the cross-section data of Ref. [21]. Indeed, an inspection 

288 of the cross-section data shown in (b) of Fig. 7 from all the works mentioned in Table III 

289 reveals that the reported cross sections are remarkably similar: the maximum difference 
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290 between cross sections is less than a factor of two; cross sections obtained separately from 

291 real and virtual photon experiments only vary by 20% to 40%. It thus appears that different 

292 methods of data interpretation, rather than cross-section determinations, give rise to the 

293 difference in reported resonance parameters. The main difference between the analyses of 

294 Refs. [21, 26] and the present analysis involves the use by the former of Siegert's theorem 

295 for extracting the B(El)i. 



296 IV. a (an, 7) ^Be RATE CALCULATION 

297 A. Reverse reaction cross section 



298 The ^Be (7, n) ^Be cross section is transformed into the ^Be (n, 7) ^Be cross-section using 

299 the reciprocity theorem. Defining ai to be the cross section for ^Be + n -> ^Be + 7, and (72 

300 to be the cross section for ^Be + 7 -> ^Be + n, the reciprocity theorem gives 

2(2jOBe + l) .^.s 
"^=(2,SBe.l)(2j„.l)^"^' 

301 where k"^ = 2yu£'„/i"^; k'^ is defined in Eq. 9; and the ground state spins for ^Be, ^Be, and a 

302 neutron are 3/2, 0, and 1/2, respectively. 

303 For several years prior to 1999, the a (an, 7) ^Be rate used in reaction network codes was 

304 adopted from Ref. [43], which considered resonant-only decays of ^Be+7 *Be + n. In other 

305 words, when considering the ^Be + n ^Be + 7 direction for the reaction, the width of the 

306 ground state of '^Be was disregarded. The rate published by Ref. [42] (known as NACRE) 

307 included the off-resonant contributions to the Q;(Q;n, 7)^Be rate. Other rates [26, 40] have 

308 since followed the formalism developed by NACRE. 



309 B. Rate calculation 

310 The derivation of astrophysical reaction rates has been described in detail in Ref. [38]. 

311 The rate per particle pair (av) is given by 

1 3 00 



312 where /i is the particle pair reduced mass, k is Boltzmann's constant, and T is the tem- 

313 perature. Equation 19 is the appropriate form for a two-body reaction. However, here the 
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-E 
kf 



EdE, 



(19) 



314 rate has to be computed for two sequential reactions. The form for calculating the rate 

315 of formation of ^Be involves constructing a double integral, taking into account the rate of 

316 formation of ^Be from the a + a scattering cross sections. The formalism adopted in the 

317 present work was developed in Ref. [44] for calculating the on- and off-resonant formation 

318 of ^^C via the triple-a reaction, and was first modified in Ref. [42] to calculate the rate of 

319 formation of ^Be for the NACRE compilation. 

320 Two a-particles interact with center-of-mass (CM) energy E to form ^Be. Subsequently, 

321 the '^Be nucleus interacts with a neutron with new CM energy E' relative to E (see Fig. 2). 

322 The rate equation has the form 



t2 /„„,\'^'-'" — AT I " ' I / H'aa ^ ^ 



with 



3/2 



X 




J a„8Be(^'; E) exp{-E'/kT)E'dE'. (21) 



324 Equation 20 is evaluated numerically using the parameters from Table II and the a + 

325 a scattering cross-sections from Ref. [45] over the temperature range 0.001 GK < T < 10 

326 GK. Table IV lists low, recommended, and high values for the presently determined rates 

327 versus temperature. The low and high rates are computed by considering the systematic 

328 uncertainty in the deconvoluted ^Be (7, n) ^Be cross section and the uncertainty in the fitting 

329 parameters (see Table II). At all temperatures, the a(an, 7)^Be rate originates primarily 

330 from the 1/2+ state in ^Be (see Fig. 9). Thus, the rate uncertainty is dominated by the cross- 

331 section uncertainty for that state. Uncertainty in the rate increases at higher temperatures, 

332 where higher states in ^Be begin to contribute noticeably to the a {an, 7) ^Be reaction rate. 
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TABLE IV: Low, recommended, and high rates for the 
a(an,7)^Be reaction versus Tg [= 1 GK ) computed from 

the present parameters. 



Tg 


Low 


Recommended 


High 


Tg 


Low 


TVT / V 

Recommended 


High 


0.001 


1.15E-59 


1.20E-59 


1.26E-59 


0.14 


4.39E-08 


4.62E-08 


4.84E-08 


0.002 


9.64E-48 


l.OlE-47 


1.06E-47 


0.15 


6.49E-08 


6.82E-08 


7.14E-08 


0.003 


6.01E-42 


6.32E-42 


6.62E-42 


0.16 


9.05E-08 


9.51E-08 


9.96E-08 


0.004 


2.76E-38 


2.90E-38 


3.04E-38 


0.18 


1.54E-07 


1.62E-07 


1.70E-07 


0.005 


1.12E-35 


1.18E-35 


1.24E-35 


0.2 


2.32E-07 


2.43E-07 


2.55E-07 


0.006 


l.llE-33 


1.17E-33 


1.23E-33 


0.25 


4.49E-07 


4.72E-07 


4.95E-07 


0.007 


4.40E-32 


4.62E-32 


4.84E-32 


0.3 


6.51E-07 


6.84E-07 


7.17E-07 


0.008 


9.21E-31 


9.67E-31 


l.OlE-30 


0.35 


8.07E-07 


8.48E-07 


8.89E-07 


0.009 


1.21E-29 


1.27E-29 


1.33E-29 


0.4 


9.10E-07 


9.57E-07 


l.OOE-06 


0.01 


1.12E-28 


1.18E-28 


1.24E-28 


0.45 


9.70E-07 


1.02E-06 


1.07E-06 


0.011 


7.90E-28 


8.31E-28 


8.71E-28 


0.5 


9.94E-07 


1.04E-06 


1.10E-06 


0.012 


4.47E-27 


4.70E-27 


4.93E-27 


0.6 


9.76E-07 


1.03E-06 


1.07E-06 


0.013 


2.12E-26 


2.22E-26 


2.33E-26 


0.7 


9.10E-07 


9.56E-07 


l.OOE-06 


0.014 


8.64E-26 


9.08E-26 


9.52E-26 


0.8 


8.29E-07 


8.71E-07 


9.13E-07 


0.015 


3.12E-25 


3.28E-25 


3.43E-25 


0.9 


7.45E-07 


7.83E-07 


8.21E-07 


0.016 


l.OlE-24 


1.06E-24 


1.12E-24 


1 


6.66E-07 


7.00E-07 


7.34E-07 


0.018 


8.24E-24 


8.65E-24 


9.07E-24 


1.25 


5.02E-07 


5.28E-07 


5.53E-07 


0.02 


5.09E-23 


5.35E-23 


5.61E-23 


1.5 


3.83E-07 


4.02E-07 


4.22E-07 


0.025 


2.51E-21 


2.64E-21 


2.77E-21 


1.75 


2.98E-07 


3.13E-07 


3.28E-07 


0.03 


4.31E-19 


4.53E-19 


4.75E-19 


2 


2.37E-07 


2.49E-07 


2.61E-07 


0.04 


1.81E-15 


1.90E-15 


2.00E-15 


2.5 


1.58E-07 


1.66E-07 


1.74E-07 


0.05 


2.63E-13 


2.76E-13 


2.90E-13 


3 


1.12E-07 


1.18E-07 


1.24E-07 


0.06 


6.88E-12 


7.23E-12 


7.58E-12 


3.5 


8.43E-08 


8.88E-08 


9.33E-08 


0.07 


6.81E-11 


7.15E-11 


7.50E-11 


4 


6.59E-08 


6.95E-08 


7.32E-08 


0.08 


3.68E-10 


3.87E-10 


4.06E-10 


5 


4.41E-08 


4.67E-08 


4.94E-08 
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FIG. 9. (Color Online) Relative contributions of each resonance to the total rate as a function 
of temperature. The temperature range relevant for the r-process, 1 GK < T < 5 GK, is shaded 
(green) . 



TABLE IV: (Continued) 







N/i(aan) 








NA{aan) 




Tg 


Low 


Recommended 


High 


Tg 


Low 


Recommended 


High 


0.09 


1.34E-09 


1.40E-09 


1.47E-09 


6 


3.22E-08 


3.43E-08 


3.65E-08 


0.1 


3.67E-09 


3.86E-09 


4.05E-09 


7 


2.49E-08 


2.68E-08 


2.86E-08 


0.11 


8.27E-09 


8.69E-09 


9.11E-09 


8 


2.02E-08 


2.18E-08 


2.35E-08 


0.12 


1.60E-08 


1.68E-08 


1.77E-08 


9 


1.69E-08 


1.84E-08 


1.99E-08 


0.13 


2.77E-08 


2.92E-08 


3.06E-08 


10 


1.46E-08 


1.59E-08 


1.74E-08 



333 Figure 10 displays the comparison of four existing rates both by their absolute magni- 

334 tudes and by normalization to the NACRE rate [42]. The present rate is 40% larger than 

335 the NACRE result for the energy range 1 GK < T < 5 GK, most important for r-process 

336 nucleosynthesis. The largest difference between the various rates exists in the off-resonant 

337 region: there the present rate is smaller than the NACRE rate by a factor of 4 because the 

338 energy dependence of all resonant cross-section contributions near the two-body threshold 

339 has been included. 
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FIG. 10. (Color Online) (a) The reaction rate as a function of temperature. The present rate is 
40% larger than the rates of Refs. [40] and [42] at the peak near T - 0.5 GK. The temperature 
range relevant for the r-process, 1 GK < T < 5 GK, is shaded (green), (b) The ratio of the rates to 
the NACRE rate. Bands shown indicate ±\a uncertainties for NACRE and the present rate. 



340 C. Rate comparisons 

341 In principle, the precision of the measured ^Be (7, n) ^Be reaction cross section should 

342 extend to the deduced astrophysical a(a?T,,7)^Be rate. The quantities used for calculating 

343 the a + a ^ ®Be rate derive from well known a + a scattering data [45], while accurate 

344 penetration factors are obtained from computed Coulomb wave functions. Recalling Fig. 

345 10, the present rate, when compared to the NACRE rate, is a factor of 3 lower at the lowest 

346 calculated temperatures, while it is 20% to 40% larger at astrophysical temperatures of 

347 interest for the r-process (1 GK < T < 5 GK). The change of the low-temperature rate is 

348 a direct result of including a realistic energy dependence for all neutron partial widths, a 

349 procedure which was not employed in Ref. [42]. The correct analytic form for s-, p-, and d- 

350 wave neutron penetration factors lowers all resonant contributions near threshold as shown 

351 in Fig. 8. At low temperatures, the energy-independent treatment leads to an artificially 
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352 inflated rate along with the non-physical dominant contributions to this rate by the j>-wave 

353 1/2" and d-wave 5/2+ states. 

354 Recently, Garrido et al. [46] have used three-body theoretical techniques to determine 

355 the direct a + a + n ^Be contribution to both the total cross section and reaction rate at 

356 energies below the two-body threshold. However, their calculated direct, three-body cross 

357 section exceeds the 93 nb experimental upper limit [37] at energies between the three- and 

358 two-body thresholds. The results presented in the present paper consider only sequential 

359 reactions and do not address the possibility of contributions by three-body processes at the 

360 lowest temperatures. 

361 The present a {an, 7) ^Be rate for 1 GK < T < 5 GK is consistently 20% to 40% larger 

362 than the rates of Refs. [40, 42]. Agreement with the NACRE rate marginally improves as 

363 the temperature approaches 10 GK. Figure 9 shows that, for the present evaluation, the 

364 1/2+ state is indeed the primary contributor to the rate at temperatures below 10 GK. Con- 

365 tributions from other states start to become noticeable for T > 5 GK. This implies that the 

366 choice of branching ratios is not important for an accurate a {an, 7) ^Be rate determination. 

367 V. CONCLUSIONS 

368 The improved accuracy of the measurement of the ^Be (7, n) ^Be reaction cross section 

369 reported in this work was made possible: a) by use of a highly-efficient neutron detector 

370 with two concentric, circular arrays of ^He tubes which provided information about the 

371 energy distribution of the detected neutrons; b) by calibration of the neutron detector's 

372 efficiency using interspersed measurement of the well-known 2H(7,n) cross section; and c) 

373 by measuring at each energy both the flux and energy distribution of the incident 7-ray 

374 beam used. Knowledge of the incident beam's energy distribution allowed deconvolution 

375 and determination of the cross section near the neutron threshold at 1665 keV. These new 

376 measurements have been used to calculate the astrophysical a {an, 7) ^Be reaction rate. 

377 Taking into account energy dependence of all neutron and 7-ray partial widths near 

378 threshold gives rise to smaller rates than previously calculated for T < 0.025 GK. For T > 2 

379 GK, contributions to the rate from higher-lying resonances become noticeable. Our cross 

380 sections and the resulting astrophysical reaction rates in the temperature range 1 GK < T < 5 

381 GK are 20% to 40% larger than previously reported. The present rate is computed using a 
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382 cross section known to ±10% at the 95% confidence level. This new rate should be employed 

383 in reaction network codes, especially those used to investigate r-process nucleosynthesis sites. 
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